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ABSTRACT 

We have used a doppler tomographic analysis to conduct a deep search for the starlight 
reflected from the planetary companion to HD 75289. In 4 nights on VLT2/UVES in 
January 2003, we obtained 684 high resolution echelle spectra with a total integra- 
tion time of 26 hours. We establish an upper limit on the planet's geometric albedo 
p < 0.12 (to the 99.9% significance level) at the most probable orbital inclination 
i ~ 60°, assuming a grey albedo, a Venus-like phase function and a planetary radius 
Rp = l.QRj U p. We are able to rule out some combinations of the predicted planetary 
radius and atmospheric albedo models with high, reflective cloud decks. 
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1 INTRODUCTION 

The existence of a close planetary companion to the GOV 
star HD75289 was first reported by Udry et al. (2000). 
. In common with 19 of the 117 extrasolar planets known, 
HD75289b is found to orbit within 0.1 AU of its parent star. 
Confirmation of the gas-giant nature of the close-orbiting 
. exoplanet, HD209458b (Charbonneau et al. 2000; Henry 
' et al. 2000), suggests that similar objects may reflect enough 
starlight to allow a direct planetary detection. Spectral mod- 
els of these "Pegasi planets" (Barman et al. 2002; Sudarsky, 
Burrows & Hubeny 2003) show that scattered starlight dom- 
inates over thermal emission at optical wavelengths. Su- 
darsky, Burrows & Hubeny (2003) suggest the high effective 
temperature and relatively low surface gravity of HD75289b 
may favour the formation of relatively bright, high-altitude 
silicate cloud decks, which act to scatter starlight back into 
space before being absorbed by alkali metals in the deeper 
atmosphere. 

Here we report the results of observations, conducted 
over 4 nights in January 2003 on VLT2/UVES echelle spec- 
trograph, aimed at detecting the starlight reflected from 
HD75289b. In Section 2 we summarise the basic physics 
which underpins our analysis, whilst sections 3 and 4 review 
the acquisition, reduction and processing of the raw echelle 
spectra, using techniques described more comprehensively 
in Collier Cameron et al. (2002) and Leigh et al. (2003). Fi- 
nally, in section 5 we use our results to place upper limits 
on the grey geometric albedo of the planet. 



2 SYSTEM PARAMETERS 

HD75289 (HIP43177) is a GO main sequence star with pa- 
rameters listed in Table 1. High- precision radial- velocity 
(RV) measurements obtained between Nov 1998 and Oct 
1999 (Udry et al. 2000) were used to identify a planetary 
companion HD75289b, whose properties (as determined di- 
rectly from RV studies or inferred using the estimated stel- 
lar parameters) are also summarised in Table 1. Since the 
planet publication (Udry et al. 2000), the orbital solution 
has been updated with new Coralie RV measurements using 
a weighted cross-correlation scheme (Pepe et al. 2002) with 
an appropriate template. 

As the planet orbits its host star, some of the starlight 
incident upon its surface is reflected towards us, producing a 
potentially detectable signature within the observed spectra 
of the star. This signature takes the form of faint copies 
of the stellar absorption lines, Doppler shifted due to the 
planet's orbital motion and greatly reduced in intensity due 
to the small fraction of starlight the planet intercepts and 
reflects back into space. 

With our knowledge of the stellar mass and the planet's 
orbital period we can estimate the orbital velocity of the 
planet V p (see Collier Cameron et al. 2002). The apparent 
RV amplitude K p of the reflected light is given by 

K p = V p sin i = 147 sin i km s _1 , (1) 

where the orbital inclination i is, according to the usual con- 
vention, the angle between the orbital angular momentum 
vector and our line of sight. 
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Star : 




Spectral Type 


GOV 1 


m v 


6.35 (0.013) 2 - 3 


Distance (pc) 


28.9 (0.46) 2 




i!iQc f'/i rA 4 
OloO \Q\J) 


M,(M & ) 


1.22 (0.05) 5 ' 6 ' 7 


R*(Rq) 


1.25 (0.05) 5 - 6 ' 7 


[Fe/H] 


0.27 (0.06) 4 


Prot (d) 


16.0 (3.0) 8 - 9 


vsini(km s -1 ) 


3.8 (0.6) 10 


Age (Gyr) 


5.6(1.0) 8 - n 


Planet : 




Orbital Period P or b (d) 


3.5091 (0.0001) 10 


Transit Epoch T (JD) 


52651.1566 (0.1251) 


K,{m s- 1 ) 


53.4 (0.6) 10 


a (AU) 


0.0483 (0.0020) 10 


M p sin i(M Jup ) 


0.455 (0.028) 10 



Table 1. System parameters (and uncertainties) for HD75289 
and its planetary companion. 1 Gratton, Focardi & Bandiera 
(1989), 2 ESA. (1997), 3 Gray, Napier & Winkler (2001), 4 San- 
tos, Israelian & Mayor (2001), 5 Fuhrmann, Pfciffcr & Bcrnkopf 
(1998), 6 Gonzalez & Laws (2000), 7 Gonzalez et al. (2001), 8 Udry 
ct al. (2000), 9 Noycs ct al. (1984), 10 Udry ct al. (2000) revised 
for this paper, 11 Donahue (1993). 

For all but the lowest inclinations, the orbital velocity 
amplitude of the planet is much greater than the broad- 
ened widths of the stellar absorption lines. Hence, lines in 
the reflected-light spectrum of the planet should be Doppler 
shifted well clear of their stellar counterparts, allowing a 
clean spectral separation for most of the orbit. 

By isolating the reflected planetary signature, we are in 
effect observing the planet-star flux ratio (e) as a function 
of orbital phase (ci) and wavelength (A), where 

e (a,A)^M^= P (A) 3 (a,A)^. (2) 

The phase function g(a, A) describes the variation in the 
planet-star flux ratio with the phase angle a, here a is the 
angle subtended at the planet by the star and the observer, 
and varies as cos a — — sin i cos <f>. As the shape of g(a, A) is 
unknown, current practice is to adopt a specific phase func- 
tion in order to express the results in terms of the planet-star 
flux ratio that would be seen at phase angle zero 

£0 (A)=p(A)^ ) (3) 

where p(A) is the wavelength-dependent geometric albedo 
and the orbital distance a is constrained by Kepler's third 
law. For reasons described in Collier Cameron et al. (2002) 
and Leigh et al. (2003), the g(a, A) we adopt is a polynomial 
approximation to the empirically determined phase function 
for Venus (Hilton 1992). 

In the event of a planet detection, analysis of the data 
should allow us to determine two fundamental properties: 

(i) K p , the planet's projected orbital velocity, from 
which we obtain the system inclination and planet mass, 
since M p sini is known from the star's Doppler wobble, 

(ii) eo, the maximum flux ratio observed where, by mak- 
ing theoretical assumptions about the nature of the planet, 
we can invoke equation (3) to constrain R p or p. 



UTC start 


Phase 


UTC End 


Phase 


^exp 


Jan 14 01:00:04 


0.396 


Jan 14 09:29:44 


0.497 


188 


Jan 15 00:56:42 


0.680 


Jan 15 09:21:32 


0.780 


173 


Jan 21 02:40:18 


0.410 


Jan 21 09:34:54 


0.492 


183 


Jan 22 03:35:30 


0.706 


Jan 22 09:41:01 


0.778 


140 



Table 2. Journal of observations. The UTC mid-times and orbital 
phases arc shown for the first and last spectral exposures secured 
on each night of observation in January 2003. The number of 
exposures is given in the final column. 

3 OBSERVATIONS 

We observed HD75289 during Jan 2003 using the blue arm 
of the UV- Visual Echelle Spectrograph (UVES), mounted 
on the 8.2-m VLT/UT2 (Kueyen) telescope at the Paranal 
Observatory in Chile. The detector was a single EEV CCD- 
44 array windowed to an image format of 2048 x 3000 15.0 
fim pixels, and centred at 475.8 nm in order 102 of the 41.6 g 
mm echelle grating, giving good wavelength coverage from 
406.5 nm to 522.0 nm. For reasons described in greater detail 
in Collier Cameron et al. (2002), we believe the prominence 
of Rayleigh scattering, combined with the high density of 
absorption lines within this spectral range, offer the best 
chance of detection. With an average pixel spacing close to 
1.5 km s _1 , the FWHM intensity in the central orders of 
the thorium-argon calibration spectra was found to be 4.5 
pixels, suggesting an effective resolving power R — 43000. 

Table 2 details the 4 nights of observations that con- 
tribute to this analysis. Each of the stellar spectra were ex- 
posed between 100s and 300s, depending upon the seeing, 
in order to expose the CCD to ~ 40000 ADU pix -1 in the 
brightest parts of the image. A typical exposure in average 
seeing (0.8 arcsec) yielded ~4x 10° electrons per pixel step 
in wavelength in the brightest orders after extraction. The 
CCD's 7s fast- readout time allowed an observing efficiency 
> 90%, and saw a gain of 2.4 and readout noise of ~ 6 e~. 



4 DATA PROCESSING 

One-dimensional spectra were extracted from the raw im- 
ages using an automated pipeline reduction system built 
around the Starlink ECHOMOP and FIGARO packages. 

After initial tracing of the echelle orders, an extraction 
process subtracted the bias from each frame, cropped it, 
determined the form and location of the stellar profile on 
each image, subtracted a linear fit to the scattered-light + 
sky background across the spatial profile, and performed an 
optimal (profile and inverse variance-weighted) extraction 
of orders across the full spatial extent of the object +sky re- 
gion. In all, 25 orders (Nos 90-114) were extracted from each 
image, giving good spectral coverage from 406.5 to 522.0 
nm. Following extraction, the S:N in the continuum of the 
brightest orders was typically 650 per pixel step. 

Flat-field balance factors were applied prior to extrac- 
tion by summing the 50-100 flat- fields taken at the start 
and end of each night. Although the UVES instrument was 
very stable between observations, we found that the noise 
was slightly reduced (~ 2%) by the use of nightly flat-fields 
rather than master flat-fields combining all 4 nights' data. 
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Simulated Planet Data Original Data 

aa oo 
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— 1 — 1 

Velocity km s Velocity km s 

Figure 1. Time series of deconvolved profiles derived from the VLT/UVES spectra, secured over 4 nights in Jan 2003. The left plot 
(a) demonstrates the result of adding a simulated planet signal at an inclination of 80° , assuming grey geometric albedo p = 0.4 and a 
radius, 1.6 Rjup- The planetary signature appears as a dark sinusoidal feature crossing from right to left as phase increases and centred 
on the superior conjunction at phase 0.5. The right plot (b) shows the original data, where no such trail is evident. The greyscale runs 
from black at — 10~ 4 times the mean stellar continuum level, to white at +10 -4 . The velocity scale is in the reference frame of the star. 



4.1 Extracting the planet signal 

The first step in extracting the reflected component is to 
subtract the direct stellar component from the observed 
spectrum. The planet signal should then consist of faint 
Doppler-shifted copies of each of the stellar absorption lines, 
at this stage buried deeply in the noise. A detailed descrip- 
tion of this procedure is given in Collier Cameron et al. 
(2002) Appendix A. After cleaning any correlated fixed- 
pattern noise remaining in the difference spectra (see Col- 
lier Cameron et al. (2002) Appendix B), we create a com- 
posite residual line profile by deconvolving each residual 
spectra with a list of stellar absorption line strengths (Col- 
lier Cameron et al. (2002) Appendix C). The deconvolution 
results in a 32-fold improvement in the S:N, by combining 
the properties of 2360 images of 1744 absorption lines listed 
within the observed wavelength range. The composite line 
profiles are then stacked by phase to display temporal vari- 
ations in brightness and RV, as at Fig. 1. Any planet signal 
present appears as a dark sinusoidal feature in velocity space 
centred on the superior conjunction at phase 0.5. 

The central pattern of distortions in Fig. 1 is consistent 
with sub-pixel shifts in the position of the stellar spectra 
with respect to the detector over the course of the night. 
Fortunately they only affect a range of velocities at which 
the planet signature would in any case be indistinguishable 
from that of the star. 



4.2 Simulated planet signatures 

We verified that a planetary signal is preserved through the 
above sequence of operations, in the presence of realistic 
noise levels, by adding a simulated planetary signal to the 
observed spectra (Fig. la). The fake signal also acts to cali- 
brate the strength of any genuine observations and the asso- 
ciated confidence contours. Detailed explanations of the sim- 
ulated planet calibration can be found in Collier Cameron 
et al. (2002) and Leigh et al. (2003). To ensure a strong 
signal we used a simulated planet of radius 1.6 Rj up and 
wavelength-independent geometric albedo p = 0.4, which 
when viewed at zero phase angle should give a planet-star 
flux ratio e = hi f* = 9-71 x 10~ 5 . 

Except for the case of a tidally locked system, any rela- 
tive motion between the surface of the planet and the surface 
of the star will result in a difference between the line widths 
of the incident and reflected spectra. From the known stel- 
lar parameters (Table 1), we have conducted Monte Carlo 
trials to estimate that the reflected component within the 
HD75289 spectra will exhibit typical line widths close to 
13 km s _1 . In order to correctly mimic the reflected starlight 
we have selected the faster rotating GOV star HD1461 as our 
planetary template. A star which is of very similar temper- 
ature and elemental abundance to HD75289. 



© 0000 RAS, MNRAS 000, 000-000 



4 C.Leigh, A. Collier Cameron, S.Udry, J-F.Donati, K.Horne, D.James and A.Penny 



Simulated Planet Data Original Data 




Figure 2. Relative probability map of model parameters K p and logeo = logp(-R p /a) 2 , derived from the VLT/UVES observations of 
HD75289. The greyscale denotes the probability relative to the best-fit model, increasing from for white to 1 for black. The contours 
show the confidence levels at which spurious detections due to non-Gaussian noise can be ruled out. From top to bottom, they show the 
99.9%, 99.0%, 95.4% and 68.3% confidence limits on eo f° r fixed K p . The left plot (a) demonstrates how the simulated planet signal is 
recovered through the analysis with ease, and is detected well above the 99.9% confidence limit. The right plot (b) containing the original 
data, however, demonstrates little evidence for a planetary detection. It is important to note that the contours only give the probability 
of a false detection if the value of K p is known in advance, which is not the case here. The slightly increased probability density seen 
close to K p = 110 km s _1 (within the K p distribution predicted by Leigh, Collier Cameron & Guillot 2003) returns a large false-alarm 
probability of 32% when considering all plausible values 35 < K p < 147 km s _1 , which is far too uncertain to claim as genuine. 



4.3 Matched Filter Analysis 

The next step involved a matched-filter analysis, described 
by Collier Cameron et al. (2002) Appendix D, to search 
for features in the time-series of composite residual profiles 
whose temporal variations in brightness and RV resemble 
those of the expected reflected-light signature. 

The relative probabilities of the \ 2 fits to the data for 
different values of the free parameters R p /a and K p are con- 
veniently displayed in greyscale form as a function of K p and 
logeo = \ogp(R p /a) . Fig. 2 details the probability map for 
the simulated (a) and original (b) observations, normalised 
to the Av 2 of the best-fit model. 

To set an upper limit on the strength of any planet sig- 
nal, or to assess the likelihood that any signal detection is 
spurious, we need to compute the probability of obtaining 
similar improvements in \ 2 by chance alone. We achieve this 
using a "bootstrap" procedure to construct empirical distri- 
butions for confidence testing, using the data themselves. In 
3000 trials, we randomized the order in which the 684 obser- 
vations were secured, but associated them with the original 
sequence of dates and times. Genuine signals are scrambled 
in phase, but re-ordered data are as capable as the original 
data of producing spurious detections through chance align- 
ments of systematic errors along a sinusoidal path through 
the data. We record the least-squares estimates of log eo and 
the associated \ 2 as functions of K p in each trial. 



The percentage points of the resulting bootstrap distri- 
bution are shown as contours in Fig. 2. From bottom to top, 
the contours represent the 68.4%, 95.4%, 99.0% and 99.9% 
bootstrap upper limits on the strength of the planet signal. 
Thus the 99.9% contour represents the value of log eo that 
was only exceeded in 3 of the 3000 trials at each K p . 



5 RESULTS 

The results of this analysis appear in the form of a rela- 
tive probability map of model parameters K p and log eo = 
logp(R p /a) 2 , shown at Fig. 2. The calibrated confidence 
contours allow us to apply constraints to the geometric 
albedo of the planet, given certain theoretical assumptions 
about the system. 

5.1 Upper Limits on Grey Albedo 

The grey albedo model makes the unlikely assumption that 
the planet-star flux ratio remains independent of wave- 
length. By adopting a theoretical radius for the planet we 
can use equation (3) to constrain its geometric albedo. Fig. 3 
details the upper limits on the albedo at the 99.9% level of 
significance, for a range of possible inclinations and radii. 

We find that the albedo limits are less strongly con- 
strained with low inclination systems. This is a natural con- 
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Upper Limits pn Grey Geemetric Alpedo - 99.9% confidence 




1 ' ' ' ' ' ' ' ' ' 1 

30 35 40 45 50 55 60 65 70 75 80 

Orbital Inclination (i°) 

Figure 3. Upper limits (99.9% confidence) on the geometric 
albedo p as a function of inclination and radius, assuming the 
atmosphere of the planet imposes no wavelength dependence on 
reflected light, i.e. a grey albedo model. The contours, from top to 
bottom, represent upper limits assuming theoretical radii for the 
orbiting planet of 1.0, 1.2, 1.4 and 1.6 Rj U p, covering orbital incli- 
nations 30° < i < 80°. If the planet exhibits the radius, 1.6 Rj up , 
and inclination, i = 75°, predicted by Leigh, Collier Cameron & 
Guillot (2003), then our observations suggest we can be 99.9% 
confident of a grey geometric albedo p < 0.12. 

sequence of the matched-filter analysis, where fitting such 
models incorporate more pixels close to the noisier ripples 
in the central stellar region of the deconvolved profile. 

If we adopt the Leigh, Collier Cameron & Guillot 
(2003) theoretical radius for HD75289b of R p = 1.6 R Jup , 
we find the planet-star flux ratio, at the most probable 
K p ~ 127 km s _1 , is constrained to f p /f* < 4.18 x 10~ 5 
at the 99.9% level, i.e. a geometric albedo less than 0.12. 
This would strongly rule out atmospheric models which in- 
corporate high-altitude, reflective cloud decks such as the 
Class V model of Sudarsky, Burrows & Pinto (2000). 

6 CONCLUSION 

We have observed HD75289 over 4 nights in Jan 2003, us- 
ing the VLT2/UVES instrument, in an attempt to detect 
starlight reflected by the known close-orbiting planetary 
companion. The excellent stability of the spectrograph, com- 
bined with the reasonable observing conditions, have en- 
abled us to produce deep upper limits on the geometric 
albedo of the planet. In truth, with our knowledge of the 
system and current theoretical predictions, we would have 
expected a reasonably unambiguous detection of the planet. 
However, we find very little evidence to suggest its presence. 
Possible reasons for the discrepancy are, 

(i) The inclination of the system is such that the spec- 
tral features of the planet do not cleanly separate from the 
stellar features, i.e. maximum Doppler shift is not sufficient 
to disentangle the planets signal. However, evidence gen- 
erated (Leigh, Collier Cameron & Guillot 2003) using the 
known system parameters, suggests an inclination i 3> 20°. 

(ii) Our choice of phase function g(ce, A) is wrong. It 
may be that close-orbiting Pegasi planets are less prone to 
back-scatter incoming starlight than we see with Jupiter or 
Venus (Hovenier 1989; Seager, Whitney & Sasselov 2000). 



(iii) Although the transiting planet HD209458 has indi- 
cated the presence of a gas giant, it may be that HD75289b 
is a more compact terrestrial planet reflecting substantially 
less starlight (see Guillot et al. 1996). 

(iv) It could be that the geometric albedo is inherently 
low (cf. Class IV CEGP of Sudarsky, Burrows & Pinto 2000), 
with starlight absorbed deep in the atmosphere of a planet 
where no high-level clouds are present to reflect the incident 
light. However, at the wavelengths observed, we do expect 
a more significant contribution from Rayleigh scattering. 

Whatever the cause, these observations provide a strong 
test for developing theoretical models which aim to predict 
the atmospheric nature of these objects. This is a field in 
desperate need of continued observations on the brightest of 
these Pegasi planets, as are detailed in the work of Leigh, 
Collier Cameron & Guillot (2003). 
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